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Abstract. Final state radiation (FSR) in pion-pair production cannot be calculated reliably because of 
the composite structure of the pions. However, FSR corrections have to be taken into account for a precise 
evaluation of the hadronic contribution to p — 2 of the muon. The role of FSR in both energy scan and 
radiative return experiments is discussed. It is shown how FSR influences the pion form factor extraction 
from experimental data and, as a consequence, the evaluation of ajj^'*. In fact the 0{a) FSR corrections 
should be included to reach the precision we are aiming at. We argue that for an extraction of the desired 



FSR-inclusive cross section o-^2i ^ photon-inclusive scan measurement of the "e 



e TT+TT -I- photons" 

cross section is needed. For exclusive scan and radiative return measurements in contrast we have to rely 



on ad hoc FSR models if we want to obtain either u^J^ or the FSR-exclusive cross section u^^^- We thus 
advocate to consider seriously precise photon-inclusive energy scan measurements at present and future 
low energy e^e~-facilities. Then together with radiative return measurements from DA^NE and BABAR 
and forthcoming scan measurements at VEPP-2000 we have a good chance to substantially improve the 
evaluation of a?,^'^ in the future. 
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^ • 1 Introduction 



Photon vacuum polarization effects are sizable and there- 
fore play an important role in electroweak precision physics. 
Because of the strong interactions between quarks and glu- 
ons the contributions of the low energy hadrons cannot be 
calculated by perturbative QCD. However, they may be 
obtained via a dispersion integral over the experimental 
e'^e~ annihilation data. A precise evaluation of hadronic 
effects in quantities like the running fine structure con- 
stant a(s) and of the muon anomalous magnetic moment 
thus depends directly on the precisioiT_of k)w energy 
"e+e" 



* hadrons" cross sections (Thad [|1|7|2H^- Further 
theoretical efforts may help to some extent to reduce the 
theoretical uncertainties of these quantities [Q. However, 
new measurements of Chad are indispensable for achieving 
substantial progress. Indeed, remarkable improvements have 
been achieved in recent years by the Collaborations CMD- 
2 at Novosibirsk and BES-II |^ at Beijing. New results 

* Work supported in part by TMR, EC-Contract 
No. ERBFMRX-CT980169 (EURODA<?NE), EC-Contract 
No. HPRN-CT-2000-00149 (Physics at Colliders), TARI 
project HPRI-CT-1999-0008 and by the Alexander von 
Humboldt Stiftung. 



are expected soon from radiative return experiments by 
KLOE (7[|l at the ^--factory DA^-NE at Frascati and from 
the B-factory at SLAG with BABAR §. 

The muon g-2 experiment at Brookhaven now has reached 
the level of 0.7 ppm in precision [ p^[Tl| ] for a measure- 
ment of Ou and depending on which evaluation of ajj'*'^ is 
adopted |Q] reveals a deviation from the Standard Model 
prediction which could be as large as 3 standard devi- 
ations. Since the main source of uncertainty of the SM 
prediction arises from the hadronic contributions, a care- 
ful reconsideration of the determination of a]^,'^'^ is manda- 
tory. In fact existing low energy "e+e^ — > tt+tt^" data are 
inconsistent with the corresponding 1=1 part obtained 
via CVC (conserved iso- vector current) from hadronic r- 
decay spectra This is a problem which most likely 

can only be resolved by new experiments. Needless to 
say that the experimental inconsistencies also reduce our 
possibilities to obtain a more precise determination of 
ajj^'^ and hence to draw conclusions about possible "new 
physics" which also would contribute to a^. 

Experiments that measure Chad do this either via an en- 
ergy scan (2m^ < < -Emax) or they measure the invari- 
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ant mass distribution of the hadronic final states ddhad/c^s' 
(s' < s) at meson factories running at fixed s, using the 
radiative return due to the emission of hard initial state 
photons. From dahad/ds' the cross section crhad(s') can 
here be extracted by factoring out the photon radiatiorj^. 

At increasing precision it becomes more and more difficult 
and challenging to extract the relevant "pseudo observ- 
able" quantities with adequate precision. By "pseudo ob- 
servable" one understands quantities obtained from raw 
experimental data only via some theoretical input. For 
example, one has to unfold the raw data from photon 
radiation effects, where the initial state radiation (ISR) 
is universal to all e"'"e~-annihilation processes, while the 
final state radiation (FSR) and the initial-final state in- 
terference (IFS) are process specific. The "pseudo observ- 
able" we are interested in is the "hadronic blob" which 
corresponds to the imaginary part of the correlator of two 
hadronic electromagnetic currents: the one-photon irre- 
ducible contributions to the photon vacuum polarization 
(see and references therein). Here we concentrate on 
low energy 7r-pair production, a relatively simple hadron 
production channel which is dominating the hadronic con- 
tribution to the muon g-2. 

For both the scan and the radiative return method we 
are facing three major sources of uncertainty affecting the 
extraction of (7^^: the experimental error, the theoretical 
error due to neglecting higher order QED corrections and 
finally the uncertainty related to non-perturbative effects 
related to photon radiation from the final hadronic state. 
Currently, great efforts are made to reduce the experimen- 
tal error below the one per cent level |||,@,||,||. QED cor- 
rections concerning low energy pion pair production have 
been considered e.g. in [^|lj, 1^,^. In the present article 
we will focus on the last of the mentioned error sources, 
the model error related to photons radiated from the final 
hadronic state. Here the problem is that the radiation of 
photons by the pions is poorly understood theoretically. 
Since perturbative QCD breaks down at low energies it 
is not possible to treat the final state pions in terms of 
their constituent quarks. On the other hand hard photons 
participating in the scattering process do probe the pion 
sub-structure. Treating pions as point-like scalar parti- 
cles by simply applying scalar QED (sQED) is therefore 
also not a solution to the problem. What makes things 
even more complicated is the fact that we have to deal 
with non-perturbative QCD effects like intermediate p or 
to resonances and photon radiation from such a hadronic 



^ As has been pointed out in already, the radiative re- 
turn "mechanism" at leading order has the nice property that 
the usual convolution integral, relating the observed cross sec- 
tion (which includes photon radiation effects) to the physical 
cross section of actual interest, appears de-convoluted (photon 
radiation acts as a spectral analyzer) such that instead of fac- 
torization under convolution integrals one has point by point 
factorization. Higher order effects which give raise to multiple 
convolution integrals of course spoil this simple picture since 
by taking one derivative we get ride of one integration only. 



state cannot be treated in a straightforward way. How- 
ever, this contribution of real photon emission can be ex- 
pected to be less important than the radiation from the 
final state pions. This is because the net charge of the in- 
termediate hadronic state is zero and in addition the de 
Brogiie wavelength of the dominant p and to resonances 
is relatively small in respect to the typical wavelength of 
the radiated photons. Only sufhciently hard photons are 
able to probe the sub-structure of a hadronic composite 
state. The pions on the other hand are charged and have a 
much longer de Brogiie wavelength. Unfortunately a simi- 
lar argument does not help for the virtual corrections since 
virtual hard photons are always included and also cannot 
be eliminated by cutting out the "trouble-making" part 
of the phase space as it is possible for real photons. As a 
consequence their magnitude is not known and they can- 
not be subtracted from the hadronic final state without 
relying on specific models like sQED. 

As mentioned above the quantity of interest is the correla- 
tor of the hadronic component of two electromagnetic cur- 
rents including strong as well as electroweak corrections. 
From the latter only the photonic corrections are sizable. 
They correspond to the FSR correction in the hadron pro- 
duction processes ||T^,Q. Since these corrections cannot 
be calculated reliably, the aim is to measure, if possible, 
the hadronic cross section a^^J^ that is dressed by final 
state photons. We would obtain in this way directly the 
quantity to be inserted into the dispersion integrals for the 
hadronic contribution to the running fine structure con- 
stant a{s) and to the muon anomalous magnetic moment 
a^, respectively, at the next to leading level of accuracy. 

Including FSR means to include photonic corrections to 
the irreducible hadronic photon self-energy. We thus ad- 
dress the question whether we can circumvent the FSR 
problem by performing an inclusive measurement, i.e., un- 
dress the data from ISR only. The question has been dis- 
cussed already in a previous paper for the radiative return 
scenario [ p^ . The result was that in this case an inclusive 
measurement does not yield the quantity of interest at 
sufficient precision. In other words, without substantial 
loss of precision one cannot avoid the necessity to undress 
from all photon radiation, including the complete treat- 
ment of FSR. Lacking a precise theoretical understand- 
ing of photon radiation by hadrons, however, one has to 
rely on model assumptions like "generalized sQED" for 
the pions (treating pions as point-like modulo a form fac- 
tor) to extract the undressed (FSR-exclusive) cross section 
CT^^^^(s) from the data in a first step. In order to obtain 

the FSR-inclusive cross section crhad(*) one has to add the 
appropriate FSR contribution "by hand" at the end. 

On the other hand for completely inclusive scan measure- 
ments we can use the fact that ISR and FSR factorize 
to "subtract" ISR from the observed inclusive total cross 
section dobs, leaving, up to 0{a^) IFS contributions, the 
desired FSR-inclusive cross section cr^2ii^)- ^^^^ 
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see for pion pair production such a "subtraction" of ISR 
is also possible with excellent precision for realistic cuts on 
the pion angles provided they are chosen such that they 
break the ISR(E)FSR-factorization only slightly. The in- 
clusive measurement requires a high quality detector with 
high acceptance and good separation of 7r° vs. 7 (Tr+Tr^Tr" 
background) . 

At CMD-2 H so far a different strategy has been used 
which we will call the exclusive scan measurement. Here 
an event selection is applied such that only soft real pho- 
tons are included which then can be corrected away. While 
there are no problems with real hard photons in this case 
one still has the problem that the virtual photon contri- 
butions from the loops include hard photon effects. The 
virtual contributions must be subtracted as well as one 
has to apply the Bloch-Nordsieck construction in order to 
get an infrared-finite cross section. Because the effective 
theory applied (generalized sQED) is renormalizable one 
obtains infrared (IR) and ultraviolet (UV) finite results. 

Note that existing data at present do not allow us to de- 
termine o'^ad('*) required for a precise determination 
of its contribution to a^'^'^. Modeling FSR by sQED we 
may estimate the size of the effect we have in mind: it is 
given at leading perturbative order by 6^a}^^ = (38.6 ± 
1.0) X 10~^^. This has to be confronted with the final preci- 
sion 6°^^a^^'^ ^ 40 X 10~^^ expected from the Brookhaven 
muon g-2 experiment. 

Before settling this 1 a (in terms of the expected final 
experimental precision) effect, it is urgent to clarify the 
origin of the 3 times larger discrepancy between e~^e~ — > 
7r+7r~ data and corresponding data obtained via CVC 
from T spectral-functions (in the energy range just above 
the p-resonance) , and the present unclear status of the 
p mass and width ||l9t] . This issue can certainly be set- 
tled by the radiative return experiments with KLOE |^ at 
LNF/Frascati and with BABAR g at SLAC0. However, 
a new energy scan experiment, which is anyway manda- 
tory for a clean measurement of the FSR-inclusive cross 
section, could also help to clarify the origin of the observed 
deviations. In addition, as we shall argue below, in radia- 
tive return measurements in order to get rid of the model 
dependent FSR contribution, at least one of the follow- 
ing conditions has to be fulfilled: i) (T7r-n-(s' < s) 3> cr^7r(s) 
(true especially for the p resonance region); ii) s' ~ s (soft 
photon region); iii) suppression of FSR by kinematic cuts. 
As we will see, at (/)-factories model dependence becomes 
an insurmountable problem at low energies below about 

^ We should mention that another very problematic energy 
region exists where experimental data are very poor or even 
controversial and which is important for the precise evaluation 
of a'^'^: the energy range 1.4 to 2.0 GeV (between the upper 
limit of the VEPP-2M machine at Novosibirsk and the lower 
limit of the BEPC machine at Beijing). Radiative return mea- 
surements with BABAR and results expected from VEPP-2000 
(upgraded VEPP-2M) will substantially improve results in this 
range. 



500 MeV where we have to deal with large contributions 
of hard FSR photons which cannot be suppressed by cuts. 

The above remarks together with the results presented 
in this paper strongly suggest that it would be desirable 
to revitalize the idea to perform an energy scan at the 
DA<?NE machine at Frascati in a second step after 
running as a ^-factory. 

In the next section we discuss the model error of pion 
form factor extraction connected to the radiation of pho- 
tons from hadronic final states in inclusive and exclu- 
sive scan scenarios. In addition we analyze to what ex- 
tent kinematic cuts on the pion angles change the picture 
and, for the inclusive scan scenario, consider the impact 
of the model uncertainty on the determination of a^'^'^. 
Section 3 is devoted to the model uncertainty of extract- 
ing the pion form factor [.^^(s')!^ for fixed s in radiative 
return experiments. We also address the question if the 
FSR contribution and its related model error can be esti- 
mated from a measurement of the pion forward-backward 
asymmetry Aps- In Appendix A we derive a general, 
model independent formula for the inclusive cross section 
(T(e+e~ — > 7* — > X -I- photons), where ISR and FSR are 
treated in a factorized form, X being an arbitrary non- 
photonic final state. In Appendix B some of the used for- 
mulas connected to FSR within sQED or fermionic QED 
(fQED) are collected. 



2 Model errors for inclusive and exclusive 
measurements in scan experiments 

2.1 Inclusive scenario 

We first present a case study of "e+e~ tt+tt" -|- n 7" 
(n = 0,1,2,.. .). Experimentally on an event by event ba- 
sis it is not possible to distinguish a final state from an 
initial state photon. In an inclusive measurement events 
with any number of (initial and/or final state) photons are 
counted. The major question will be to what extent and 
at what accuracy we may evaluate FSR-inclusive cross 
sections from the experimental data. We first consider the 
measurement of the FSR-inclusive cross section^ (t''''''(s) 
in energy scan experiments. 

Suppose for the moment, that we would be able to cal- 
culate photon radiation from pions. Then we would have 
two possibilities: 

(i) determine the undressed cross section ct'*'' (s) by un- 
folding the observed cross section from all photon ra- 
diation and add the FSR as calculated by perturbation 
theory with desired accuracy, which yields a^'^^; 

In the following we drop the labels "had" or "tttv" for cross 
sections like cr'"' and a'^^'K 
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(ii) determine an FSR-inclusive cross section by unfold- 
ing only the calculated ISR from the observed cross 
section, which yields a^'^^s). 

The question then is to what accuracy does a^'^^{s) ap- 
proximate a^''^{s). Since, actually, we do not know how to 
calculate photon radiation from pions in a model indepen- 
dent way only the second approach is able to give a model- 
independent answer, however, then we do not know how 
well ct'-''''' (s) approximates the quantity of interest (T^'''-'(s). 
What we will do in this case is a "guesstimate" of the qual- 
ity of the approximation by modeling FSR by generalized 
sQED. 

For the error estimate the following factorization theorem 
is crucial: Neglecting the IFS contribution, being of 0{a^) 
due to charge conjugation invariance, the inclusive total 
cross section crobs(s) may be written in a factorized form, 

CTobs(s) = J dsv (J^'^Hsv) P^nfis, sv) + 0{a%FS , 

(2.1) 

which means that FSR and ISR can be treated indepen- 
dent ly from each other. Details are given in Appendix A 
[see ( A.2lJ), (JA.28 )]. It is worth to stress that the powerful 
identity ( |2.lD is not easily recognizable to fixed pertur- 
bative order. Note that ( p.l[ ) is quite general once IFS is 
neglected. It is based on the fact that we have a neu- 
tral current process for which wc can apply a separa- 
tion into Lorentz-covariant and individually gauge invari- 
ant initial and final state contributions. Qualitatively the 
result may be understood as follows: by the fact that 
at low energies the single virtual photon exchange (1/s- 
enhancement) highly dominates the "e+e~ 7r"'"7r~ -|-n7" 
cross section and due to the suppression of the IFS (see 
later) it makes sense to consider the process in an approx- 
imation of an s-channel single photon exchange (i.e. dia- 
grams which factorize into two disconnected parts upon 
cutting the photon line) . This virtual photon then carries 
the invariant mass ^/sv and the above convolution is ex- 
act up to the indicated missing IFS effects. We would like 
to remind the reader that the representation of the tttt- 
production cross section in terms of the pion form factom 



.(0) 



(.) = |F(°)(.)PctO-p°'"*(.) 



(2.2) 



with [Pt^ = (1 — 4m^/s)-'^" is the pion velocity] 

^0, point _ TT Q^/?^ 

3 s 



* Note that (j'-'''(s) and equivalently |-F7r"''(s)p are not mea- 
surable quantities, as we shall discuss below. They are useful, 
theoretically motivated concepts defined in a world where the 
electroweak interactions are switched off. In reality we cannot 
switch off QED effects and this is part of the problem we are 
dealing with in this paper. If one could calculate l-Fi"' (s)|^ 
for time-like s non-perturbatively in lattice QCD, this is the 
quantity what one would take from lattice QCD. 



also makes sense strictly only for the one-photon exchange 
approximation]^. In this approximation sy can be neatly 
identified with s in |_F!n.(s)p, in spite of the fact that sy, 
as the squared invariant ma ss of a virtual state, is not 
an observable. Thus in ( |2.lD sy is only a formal (un- 
physical) integration variable where the boundaries are 
physical observables: 4toJ < sy < s. Nevertheless, we 
can always fit the pseudo observable a^''\sv) to the ob- 
served data crobs(s) by using (2T) and thereby determine 
(7^''\sv)- The accuracy with which this can be achieved, 
up to IFS contributions, is limited by our knowledge of the 
initial state radiator function p\^f{s,sv) only. The latter 
can be calculated without any model dependence within 
perturbative QED (see e.g. |^l|). Since IFS effects are of 
0{a^), the model dependence for the extraction of the 
FSR-inclusive cross section is determined by an (as yet 
unknown) 0{a^) IFS contribution]^. What is very impor- 
tant is that IFS interference does not include contributions 
from leading logarithms of the kind log(s/mg). We may 
estimate the 0{a^) effect to be at the per mill level. 

As we have already stressed, our "master formula" ( ^t] ) 
cannot be directly applied to a real experiment with some 
cuts and/or detector inefficiencies (the leading uncertain- 
ties are due to the need of extrapolation to the blind 
zones of the measurement). In a real experiment the in- 
fluence of these effects will be taken into account using 
a realistic Monte Carlo event generator which features a 
high quality ISR matrix element and some modeling of 
FSR. Let us focus therefore on a situation where angu- 
lar cuts are present. Then, ISR and FSR phase space in- 
tegrations can not b e disentangled as it is possible with- 
out cuts [see (ATI)] and ISROFSR factorization breaks 
down at the 0{a) level. As a consequence, to subtract 
only ISR from the data we have to rely on specific FSR 
models. To be able to extract a^^^sv) without signifi- 
cant model dependence the condition that the applied cuts 
break ISR(8)FSR factorization only slightly has to be ful- 
filled. Here we will investigate the breaking of ISR^FSR 
factorization by some semi-realistic C-symmetric cuts, 
©TTi > 6**^, G^TTi being the laboratory angle between the 
pion momenta and the beam axis, treating FSR by sQED. 
For such cuts we then can write the observed cross sec- 
tion as (for details see Appendix B and ]]l^ ; A is the soft 
photon energy separating soft from hard photons) 

s-2^A 



4m?. 



dsv aitisv)pfi^ns,sv)-6lT{s) 



(2.3) 



with Siai{s,A) corresponding to the soft plus virtual and 
pfniis, sy) corresponding to hard photon initial state QED 



Since in (2.2) FiP\s) does not d epend on the pion pro- 



duction angle a similar formula (B.18) is valid for the case of 
angular cuts with the same function ^^"'(3). 

^ The O(a^) IFS is complete for the real photon emission. 
However, some virtual contributions are not yet calculated. 
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corrections. The 0{a) ISR(8)FSR factorization breaking 
term 5^^j"(s) accounts for the missing pion events which 
cannot be seen in the experiment: 



- JO) 



(s)- {v{s)~Vcut{s)} + 0{a').{2A) 



'^cut" vanishes for the case without cuts (restoration of 
factorization). Remember that for C-symmetric angular 
cuts the 0{a) IFS contribution drops out. 



In a world with point-like pions we could calculate 6'; 
perturbatively in sQED, where 



scan 
cut 



is) 



'7(cut)(s) 



(5fi„(s,yl) + 



ds' ptT\s,s') 



4mi 



(2.5) 

with s' being the square of the invariant mass of the pion 
pair and 5f^-n{s,A) and pgj"'''(s, s') the corresponding FSR 
corrections given in Appendix B. For real world pions we 
may estimate this term assuming generalized sQED which 
at least treats the soft photon part correctly and for the 
rest is a guess. It means that we assume that (^) and 
(2.5), with ri{s) calculated in sQED, still to some approx- 
imation accounts for the effect. What we will actually do 
is to consider (5^„t"(s), evaluated as just described, as a 
theoretical uncertainty (model error). 



Note that ( |2.3| ) only contains the measured cross section 
(T^5J*(s), the known initial state correction factors Jini and 

Pi™*\ and the FSR-inclusive cross section o•^^^(s), which 
is the quantity to be extracted from the data since it cor - 
responds to the FSR-inclusive pion form factor via ( B.2Cl| ). 



Whether the approximation (T^^|(s) obtained via 



ds'ai-:^,is')p^:^is,.s'), (2.6) 



after neglecting (5^^™(s) in (pT^), yields a good approxima- 
tion for the FSR-inclusive cross section is subject of the 
investigation described in the following. 

We first introduce the FSR-inclusive form factor Fjf \s) 
by 

|Fi^)(.)P = |F(°)(.)p (l + ^^(.))+0(a2), (2.7) 
and assume that in some approximation it makes sense to 



write formulas like (2.2) also between a^^^s) and Fi''\s) 

and between a^'^'' (s) and Pi''^ (s). Hard photon effects spoil 
these assumptions at some level, but this at the moment 
is difficult to quantify. So in the following, this will be part 
of our model assumption (see below). 

To estimate the model dependence for the extraction of 
^cutis) we first generate a sample f7Q5J*(s), using the pion 



form factor IF, 



(0), 



as giv en in Q and the relations ( 2^ ) 
(|B.21| ). Then we utilize the MINUIT 



- (^ and (P.ISD 
package ||2^ to obtain 
(T^„t(s)] from the croU*(s) data using (2^). For the data 



F^''\s)\'^ [which corresponds to 



fitting we adopt again the Gounaris-Sakurai type param- 
in the version proposed in Q. 



eterization of |Fi''^(s)p 



We then estimate the model error by 



^cur(s) 



Fi-'^s) 



(2.8) 




0.6 0.8 
[GeV] 



Fig. 1. Estimated relative model error for the extraction of 
the absolute square of the FSR-inclusive pion form factor, 
|Kr^'(s)p, as a function of the center of mass energy in a 
photon-inclusive scan experiments for different C-symmetric 
angular cuts on the pion angles. Curve (aso) and curve (aeo) 
corresponds to the cuts 0^ > 30° and ©tt > 60°, respectively. 

The model error for the extraction of |Fi''''(s)p from scan 
data with C-symmetric cuts is shown in Fig. |l|. The de- 
tailed shape of the curve depends on the parameterization 
of the pion form factor. It can be noticed that the consid- 
ered cuts do not lead to large model errors. Even for the 
extreme cut of Oj^ > 60° the model error is still below half 
a per cent, for 0-^ > 30° it is below 1 per mill. In fact, the 
observed smallness of the model error is related to the p- 
wave-like angular distribution of the outgoing pions. Fig. |^ 
shows that angular cuts of the pion angle against the beam 
axis up to 30° decrease the cross section very little. 



In Fig. H the impact of the discussed model error on 
V 37r / 



^had 



^had 



ds ^^^^ , (2.9) 

4m?. * 



is shown. For this we compare the values of a^^'^ when 
inserting 

i?..(s)^i?(;)(s) = ^|F(^)(.)P (2.10) 
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a 



500 




0.6 0.8 
Vs [GeV] 



Fig. 2. Observed total cross section for pion pair production 
as a function of the center of mass energy (^/s < 1.02 GeV) 
for the considered cut scenarios. 



into (2_^) [this value is denoted by a^^'^'"''''] with the value 
when inserting 

R^s) EE ^ ^ \F(-r)(sr , (2.11) 



[this value is denoted by a^^ Then the model error 
of aj;^^ which is plotted in Fig. |[ for the curves (039) and 
(ago) is defined as 



had(7) _ ~had(7) 

/I had _ ^ 

A' ~ had (7) 

a,. 



(2.12) 



o 

X 2 



to , 



. . 1 . . . 1 . . 


1 ' ' ' 1 _ 

(360^''^- 




/' . 






- , , 1 , , , 1 , , , 1 , , , 1 - 


0.4 0.6 


0.8 1 



^/s„,„ [GeV] 



Fig. 3. Estimated relative model error in per mill of 
i)j^'*(smin < s < MJ) in scan experiments. The curves (aao) 
and (060) correspond to the cases in Fig. 0. 



Let us remind that the present theoretical error is at the 
level of 1.2% percent. 

As a first summary we may say that the direct extrac- 
tion of |Fi'^''(s)|^ in an inclusive scan yields a very good 

approximation of |i^i''^'*(s)|^, especially in the low energy 
region where the contribution of FSR becomes large. 



2.2 Exclusive scenario 



Scan measurements in the past attempted to extract the 
"bare" cross section (T^°)(s), undressed from photon radi- 
ation effects. As already mentioned, the bare cross section 
is the object of primary theoretical interest, in principle. 
It is the quantity which allows us to extract the pion form 
factor which encodes the strong interaction structure of 
the pion in a world where the electromagnetic interaction 
has been switched off. It is the non-perturbative quan- 
tity which one would compute by a simulation in lattice 
QCD or investigate by means of general low energy prop- 
erties of the strong interactions like chiral perturbation 
theory, locality and analyticity The theoretical con- 
cept of disentangling effects from different interactions has 
been very successful, however, has its limitation at some 
point. In phenomenology of low energy hadrons it is in 
fact not possible to separate in a model-independent way 
QED from QCD effects (at the level of accuracy we are 
considering here weak interaction effects are negligible). 
This will be discussed in more detail in the next section. 
Here, for the moment, we assume that (t*^°)(s) is a sensible 
pseudo observable. 

Fig. ^ shows the relative deviation of the FSR-inclusive 
pion form factor |i^i''''(s)|^ (being the desired quantity to 
be inserted into the dispersion integrals) from the un- 
dressed pio n form facto r LF^(s)|^, as calculated within 



sQED [see ( pTl ,( |BT^ ),( p!20| )] 



r(7) 



(«) 



(s) _ |Fi^)(s)i^-|Fr(s)i 



(0)/ 



cr(0)(s) 

-77(s) + 0(a2). 

TT 



Fi")(.)|2 



(2.13) 



This quantity can be taken as a measure of the importance 
of FSR. 



Since the analysis presented so far does not account for 
the fact that generalized sQED describes correctly the soft 
photon part of the FSR-spectrum onlyQ, we would like to 

^ Note that for s 4m^ real FSR is known precisely since 
there is only enough phase space for the radiation of soft pho- 
tons and for soft photons the FSR radiation mechanism is uni- 
versal for given masses and charges of the final state particles. 
Therefore estimating the model uncertainty to be given by the 
sQED result appears to be too crude as it overestimates the 
model uncertainty in the soft photon region. 
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allow us to describe also hard photons in a realistic man- 




0.4 0.6 0.8 1 < 

Vs [GeV] 



Fig. 4. Importance of FSR in sQED. The curve shows the 
difference between the absolute square of the undressed and of 
the FSR-inclusive pion form factor [see (2.13)]. 



0.4 



0.8 



y/s [GeV] 



go further and compare sQED modeling with a second 
one which differs from it for hard photons. As in JTsI we 
compare two different treatments of final state corrections: 
once we take for pfin(s,s') and (5fin(s,yl) the functions re- 
lated to photonic radiation from point-like, scalar pions 
and once we use the corresponding functions for the pho- 
tonic radiation from point-like, fermionic pions with the 
same charge and mass (see Appendix B). The reason we do 
this is that in the soft photon limit the scalar as well as the 
fermionic approach yields the same correct result for the 
real photon contribution. For hard photons on the other 
hand both scenarios are obviously different. This will al- 
low us to get some feeling in what kinematic regions hard 
photons play a substantial role and there lead to a large 
model uncertainty. 

We thus consider in the following the model dependence of 
FSR effects, by replacing the integrated final state correc- 
tions for scalar particles, represented by the factor ri{s), 
by a corresponding factor for a fermionic final state, which 
we denote by ?y'^(s) (see again Appendix B for explicit 
formulas). Though both 77(5) and ?7-^(s) diverge when ap- 
proaching the pion pair production threshold (Coulomb 
polcQ), their difference in this limit is a small number 
lims^4^2 [ti{s) — rjf (s)] = 1 (in units — )|^. As already men- 
tioned, obviously, the generahzed sQED/fQED modehng 
of FSR obtained by replacing the point-pion from-factor 
"1" by a form factor function |F7r(s)p of one single vari- 
able s is valid for soft photons only. A method which will 



The Coulomb resummation has been considered in |13 
At high energies scalar and fermionic FSR read 



7]{s 



3 and rj\s 00) = ^ 



Fig. 5. Model estimate of the relative model uncertainty for 
the extraction of the absolute square of the FSR-inclusive pion 
form factor, \F}f'\s)\'^ , as a function of the center of mass en- 
ergy in scan experiments. 



Our discussion thus motivates the consideration of the fol- 
lowing measure of the model dependence: 



/iscan 



r(7) 



a 



-[ri{s)-il^{s)] + 0{a'). 



(2.14) 



It compares in a ratio \F^""{s)\'^ , being the pion form fac- 
tor dressed by scalar FSR, with the absolute square of the 
pion form factor \F^'^^'^ {s)\'^ ^ being dressed by fermionic 
FSR and corresponds to our ignorance of FSR. Fig. || 
shows AY^^[s) as a function of energy and suggests an 
uncertainty below 0.5% over the whole energy range of 
interest. 

For the data analysis at the CMD-2 experiment |^ the 
event selection was such that only events containing real 
low energy photons were taken into account and thus the 
pions were approximately back to back. Since for soft pho- 
tons the FSR mechanism is known (factorization) , the real 
photonic corrections together with the universal soft plus 
virtual IR terms can be subtracted from the observed cross 
section in an essentially model independent way|^. Ac- 
cordingly, in order to keep the formulas simple, we define 
a subtracted cross section cr^^g^'Xa) which does not depend 
on the cuts any longer and is obtained from the experi- 
mentally observed exclusive cross section in a theoretically 
well controlled manner. 



Collinear hard photons can be easily separated via the 
event shapes. 
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While real hard photons may be eliminated by appropriate 
cuts, the same cannot be done for the remaining virtual 
corrections which include high momentum scales in loops 
and hence their treatment is model dependent. As a con- 
sequence the determination of the undressed cross section 
(t(°^ (s) as well as of the FSR-inclusive cross section a^'^^ (s) 
from the given data suffers from model dependence which 
is hard to estimate. 

Here we will present a model error estimate assuming that 
the observed cross section can be written as a product of 
(T^°-' (s) containing all QCD effects (pion form factor) and 
a function containing only the initial and final state QED 
corrections. This ad hoc assumption, although criticizable, 
seems to be the best we can do so far. Applying the proce- 
dure of real and IR photon subtraction as described above 
then yields 



^subtr / \ 



G^) l + (2.15) 



V+Sf 



which is the 0{a) cross section including only the non-IR 
initial and final state soft plus virtual corrections corre- 
sponding to S^'^^ (s) and S'^^^ (s). Here of course the final 
state correction factor S'^^^^is) is not known. Thus at this 
stage we have two unknowns: cr'^^ (s) and the FSR cor- 
rection. Only after assuming that FSR is given by sQED 
or fQED we can then extract the undressed or the FSR- 
inclusive cross section from cr^^g'''(s) via the following for- 
mulas (see Appendix B): 



(2.16) 



subtr 
obs 



(^) 



l + f 77(/)(s) 



(2.17) 



Using ( 2.16 ) we could try to estimate the uncertainty for 
the extraction of (j^"^ (s) via 



a: 



cxcl. 



i-(o) 



r(0),/ 



a(o)(i 



1 



(^) 



(2.18) 



The such estimated model error is shown in Fig. ^ [curve 
(0)]. We would like to stress that we should be careful 
not to take this error estimate obtained from the compar- 
ison of two factorizable models too seriously. In fact, if we 
would make the analogous comparison for the extraction 
of a'^'^^s) the such obtained error would be of the level of 
1 per mill. However, we would expect a larger error from 
non-factorizable FSR contributions which cannot be esti- 
mated. 

As an alternative possibility we may try, in the spirit 
of (2.1), what we get if we just correct for the model- 
independent ISR and the model-independent IR-sensitive 



part of FSR, obtaining 



is) 



-subtr 
^obs 



is) 



(2.19) 



Defining 



a: 



cxcl. 



Xs) 



is) 



is) 



a(^)(s) 
[l + ^vis)][l 



sV+S 



= 1 - 



is)] 



^ + siris) + sL^'is) 



(2.20) 



we can get a feeling for how well a^'^'^ (s) approximates the 
true a'^'^\s) [estimated here by ( 2.17 )]. The result is shown 
in Fig. H [curve (7)]. 

To summarize: what can we get from a hard-photon ex- 
clusive measurement: 

i) cr*^°^ : in spite of the fact that all real hard photons have 
been eliminated by cuts a surprisingly large model un- 
certainty due to hard virtual photons poses an inherent 
limitation: the corresponding uncertainty cannot fall be- 
low the level of about 0.5% (sQED). Strictly speaking cr(°) 
is not accessible to experiment or only at limited precision 
by the fact that we cannot switch off virtual QED effects 
in reality. 

ii) (t('^): the missing real hard photons must be calculated 
from a model like sQED and added by hand. What we 
get is a model dependent a^'''^'^'""'*"'^ (s). Surprisingly, the 
model dependence we estimate by our method (assuming 
factorization with a single scale form factor) for this ob- 
ject is much smaller (at the level of 0.1% only). On the one 
hand this reduced model dependence can be traced back to 
the Kinoshita-Lee-Nauenberg (KLN) theorem, which in- 
fers that radiative corrections for total inclusive cross sec- 
tions are free from large logs. On the other hand it is not 
conceivable that we get a more precise knowledge of a^^^ 
from not measuring hard photons than from actually mea- 
suring everything. In the latter case the uncertainty shown 
in Fig. 1^ has been estimated, which, as expected, shows an 
increasing uncertainty for increasingly strong cuts. Never- 
theless, even so we think that our method of estimating 
the model dependence underestimates the error in the ex- 
clusive case, it is a quantity which is protected by the 
KLN theorem from large effects and thus is a quantity 
which seems to be under much better control than e.g. 
the bare a^^^ . How much better is hard to quantify at this 
stage. 

iii) a^'^\s): is model independent per definition but it is 
not the quantity of actual interest, as it is not a good 
approximation to a^''^ After all the hard real photons 
are missing here and again we only can get what we are 
interested in by adding the missing piece using a model. If 
we do so we end up with a^'y)'('^°<^'^^'> (g) again, up to higher 
order terms. Then we are essentially back at ii). 
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Vs [GeV] 



Fig. 6. Model error estimations for the extraction of the ab- 
solut e squ ar e of the pion form factor in an exclusive scenario 
[see dliol)] 



3 Model errors in radiative return 
measurements 



At radiative return experiments the spectral function 
da/ds' is measured where ^/s' is the invariant mass of 
the non-photonic final sta te. Let us first have a look at 
our "master formula" ( |2.l|) which is the photon-inclusive 
cross section in form of a convolution integral in the in- 
tegration variable sy, sy being the invariant mass square 
of the hadronic final state including the FSR photons but 
excluding the ISR photons. One could think that, due to 
the factorization o f IS R and FSR already on the matrix 
element level [see ( |A.5| )], it is possible to extract a^'^\sv) 
also from a radiative return measurement at fixed s. Of 



course by rewriting (2.1) we can formally get 



dcj-u 



dsv 



rmi 



(s, sv) + 0(a^)iFS- 



(3.21) 



However, we cannot extract a^'''^{sv) from the experimen- 
tal data for the simple reason that sy is not an observable. 
This can be immediately seen already for the case of sin- 
gle photon emission where we have sy = s' if the photon 
is emitted from the initial state, but sy = s if the pho- 
ton is emitted from the final state. Since on an event level 
ISR and FSR photons cannot be distinguished, sy cannot 
be obtained from the observables s and s'. The error we 
would make when identifying s' with sy is therefore of 
leading order FSR. This is one way to see that not only 
there is no way to measure the FSR-inclusive cross section 
a^'^\s) in a radiative return measurement but, in fact, we 
have to deal with an 0(1) FSR-background leading to 
an in general significant model error for the extraction of 
even the undressed cross section (t^'^^(s). In the following 
we therefore are going to investigate the model error for 
radiative return scenarios in a separate analysis. 



Taking into account radiative corrections in the approx- 
imation where only the leading single photon radiation 
from the final state is included, the observed spectral func- 
tion can be expressed as the sum of an ISR and an FSR 
contribution since the IFS contribution drops out: 

^) = is') pini(s, s') + (s) pn^s, s'). 

(3.22) 

Again we refer to Appendix B and for the explicit ex- 
pressions. Because we are interested to measure the FSR- 
inclusive cross section it is tempting to extract the quan- 
tity 



1 



Pini{s,s') 



da 
d? 



(3.23) 



which is obtained by just subtracting the I SR pa rt and 
dropping the model-dependent last term of (|3.22| ). How- 
ever, since what we really want to get is a^^\s) we have 
to rewrite ( 3.22 ) in terms of this true 0{a) FSR-inclusive 
cross section. We easily find 



with 



.^^fl4f7(0)(s) + -,7(s')a(0)(s'). 

Pini(S,s') TT 



(3.25) 



We can see that the first term of (3.25) is of 0(1) but the 
second of 0{a), such that no cancellation up to higher 
order terms is possible! The first term may be considered 
as a correction term only in the region where we have 
a large enhancement of (t^''-'(s') by the p resonance, in 
particular, in comparison to the reference cross section 
(T^°^(s) at s = M| (or M^). In addition, the mass effects 
of photon radiation by the pions vs. the ones from the 
electron-positron system in fact lead to quite some sup- 
pression of yOfin(s,s') in comparison to pini(s,s'), both of 
which are of 0{a). We conclude that, in radiative return 
experiments, a direct model independent extraction of the 
FSR-inclusive a^'^^s') is not possible at 0{a) precision in 
the naive way just considered. 

Since we try here to discuss 0{a) corrections to (t'°^(s'), 
which in the radiative return scenario is given by 



cr 



is') 



1 



/5ini(s,s') 



da 
d? 



obs 



Pfinis,s') 
/5ini(s,s') 



is) 



(3.26) 



one obvious deficiency of our starting equation ( |3.22 ) are 
t he m issing higher order corrections. In the resolved form 
( 3.26 ) we are at the 0{1) level only and thus we are not 
able to seriously address the question of an FSR-inclusive 
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measurement. Obviously we loose one order in a in a ra- 
diative return measurement. A discussion beyond leading 
order FSR wou ld be possible only if the complete next or- 
der version of ds^H) would be available. The need for go- 
ing to higher orders also leads to more problems with the 
treatment of hard photons radiated by the pions. Of course 
the limitations are coming from the fact that virtual hard 
photon emission by the pions cannot be switched off and 
our limited theoretical knowledge of the higher order con- 
tributions have their drawback for a precise determination 
of a^^\s') or cr(^)(s'). The basic problems and limitations 
discussed above for exclusive scan measurements also ap- 
ply for the radiative r eturn method. Nevertheless, a dis- 
cussion on the basis of ( 3.22 ) addresses the major difficulty 
we encounter in the attempt to measure the FSR-dressed 
cross-section in a radiative return experiment. 



200 



150 



100 




As ( 3.24 ) tells us, the FSR corre ction ^77(5'), given pre- 
cisely by the second term of ( 3.25 ), is completely lost once 



we drop (5 ''''''(s , s') in order to get the model- independent 



quantity ( |3.23 ), which means that the latter is not a very 
meaningful quantity. Therefore, in |Q we proposed to un- 
fold the raw data from all photon radiation, by modeling 
FSR by generalized sQED. The extracted undressed cross 
section a^^\s) then suffers from model dependence. From 
the preceding discussion we know what the actual pr oblem 
is: in first place we have to control the first term (3.25) 
which is suppressed by 

Pfin(s, s')/Pim(s, s') and in regions where a'^^\s)/a^-"^s') 
is small but otherwise is of 0(1). 



A measure for the relative importance of the disturbing 
model-dependent FSR term 5^^^^ is 



(^') 



r(0) 



(^') 



(7(")(S') 



(3.27) 



which is indeed a measure for the importance of FSR as 
given by the radia tor function punis, s'). It has to be com- 
pared with ( ^.13 ) which measures the FSR in the inte- 
grated fo rm (| 2.5) for the case of a scan experiment. Ob- 



viously, ( 3.27| ) definitely does not account for the fri{s') 
term in (t'^^'(s'). We may consider it, however, as a mea- 
sure for the model-dependence of the extraction of the 
undressed a^^^s'). Thus let us point out once more, it 
would be misleading to think that d^'^\s') in any sense 
would approximate cr(^^(s') to better than the 0(1) level. 
Although it includes FSR effects, it does not include the 
ones we are looking for. 

The size of FSR effects for the radiative return scenario is 
depicted in Fig. 0. 

Clearly the effect is large both in the soft (s'<s) and 
the hard (s' <C s) photon limits. It is small only where 
it is suppressed by a large a^°\s'), i.e., around the p- 
resonance. 

In fact the raise of Sr^^^ {s' , s) in the soft photon regime 
(for s' s) just means that soft photon FSR effects be- 
come large and does not imply a large model dependence. 



Fig. 7. Importance of FSR for sQED. The curve shows the 
scenario for the extraction of 1^^"' (s')|^ in radiative return 
experiments {^/s = M^) estimated by once including and once 
excluding final state corrections [See (3.27)]. 



Hence taking ( 3.27 ) as an estimation of the model error 
would be too rough. It does not yet take into account that 
sQED describes well the soft photon regime. In order to 
get a more realistic measure for the model dependence we 
have to proceed as in the previous section. 



Given the experimental distribution {■^)^^^^ in ( 3.22 ) the 
extracted a^°'>{s') depends on the unknown FSR radiator 
Pfin(s, s'), which we again model by sQ ED. A nalogously to 
our analysis for the scan scenario [see (2.14)] we compare 
the result for a scalar vs. a fermionic radiator by looking 
at 



/ir.r 



is') = 



r(0),/ 



r(0) 



cr(0)(s' 



(3.28) 



with cr'^*'-'(s') extracted assuming sQED and a^'^'>'^{s') ex- 



tract ed a ssuming a fermionic radiator p^^{s, s') [see (B.2) 



and (B.3)] 



The result is shown in Fig. |8[ We would like to stress once 
more that (3.22) does not incorporate the ISR(8)FSR and 
IFS effects, which account for an additional 0{a) model 
error contribution. Thus, Fig. |^ when taken without this 
proviso is misleading at energies there the 0(1) term is 
kinematically suppressed to be smaller than the missing, 
presently unknown, 0{a) termsj^ which are expected to 
be at the few per mill level. 



Clearly, the "scalar versus fermionic" scenario gives, as 
expected, a small model dependence for the soft photon 
region but the model error remains large for pion pair 



Their evaluation would require a full two-loop calculation 
of the process e'^e~ tt+tt". 
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factor of about 30 larger near tttt threshold than the con- 
tribution from "e~''e~ — > 7r~''7r~ + photons" . Higher order 
photonic corrections also have to be taken into account. 
Leading log photonic 0{a^) corrections here contribute 
about 0.3 per cent to the inclusive spectral function. 

Let us now apply C-symmetric angular cuts to the calcula- 
tions concerning the r adiati ve return method. Cut objects 
have been defined in (|B.18|) - (|B.21|). 



Analogously to the case without cuts in ( 3.27 ) the im- 
portance of the FSR contribution for a C-symmetric cut 
scenario can be defined as the relative deviation of 



^(7) 
"cut 



1 



da 
lb' 



(3.29) 



obs.cut 



Fig. 8. Estimated relative model error for the extraction of 
|J^i°'(s')p in radiative return experiments (-^s = M^). The 
curve describes a scenario where the model error is estimated 
by a comparison of FSR once from scalar particles and once 
from fermions of the same charge and mass. 



from the undressed cross section a, 



(0). 



,(7) 



is')- 



JO) 



is') 



J") 



(,s') 



(3.30) 



production near threshold energies where hard photons 
are involved necessarily. 

We repeat that for radiative return experiments we can- 
not see a possibility to measure the FSR-inclusive cross 
section a^'^\s), at least not in some obvious way. We only 
are able to extract the undressed cross-section a^^\s). 
At 0{a) precision even this is only possible in a model 
dependent way, since the problems are the same as the 
ones we have addressed earlier for the exclusive scan mea- 
surements. To get a^''\s) we must add up the FSR as 
given by sQED, with the drawback that we have to live 
with the model dependence as i llustr ated by Fig. ||. Since 
the disturbing second term in ( 3.22 ) is much smaller for 
a radiative return experiment at a S-factory, it seems 
that the chances to get a model-independent determina- 
tion of (T^°^(s') there, could be good for what concerns 
the theoretical uncertainties associated with FSR. It can 
be easily checked that for such measurements at a B- 
factory indeed the model error is limited by higher or- 
der FSR effects (not considered here) since the 0(1) FSR 
contribution is essentially due to the fact that cr,r7r(s = 
M|.^g)/(T^7r(s < 1 GeV)<0.04. The observed cross sec- 
tion, on the other hand, is reduced by about two orders 
of magnitude in respect to a <?-factory measurements. 
However, for the BABAR experiment at SLAC this draw- 
back is compensated by a very high luminosity which is 
about a factor of 400 larger than at the DA<PNE col- 
lider, ft is important to note that a good control of ISR 
will be required since the gap between s' and s is much 
larger than for (^-factories. In particular the singlet ini- 
tial state pair production channel "e+e~ 7r"'"7r~e+e~" 
yields the dominant contribution to the inclusive channel 
"e+e^ TT+TT^ + anything" at B-factory energies, being 
about a factor of 3 larger around the p peak and even a 
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Vs' [GeV] 



Fig. 9. Relative FSR contribution as given by sQED obscuring 
the extraction of |i^i'^'(s')p in radiative retur n experiments 



(ys — M<p) for different angular cuts [see ( 3.3C )]. Curve (a) is 
the same as the curve in Fig. ^ Curve (b) corresponds to the 
cut scenario where only events are taken into account for which 
the laboratory angle between the pion momenta and the beam 
axis is larger than 30° and the laboratory photon angle 9^ 
is restricted to a region 7° < d-y < 20°. In a similar way curve 
(c) corresponds to 6-„ > 40° and 9-y < 25°. 



Again we may "guesstimate" a model uncertainty by re- 
placing and comparing the sQED model with the fQED 
model. In the latter case we replace Pf^{s, s') by 



1 



point, f 
^0,cut 



da 
dp 



point, f 



fin, cut 



(3.31) 
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Fig. 10. Estimated relative model error in per mill for the 
extraction of l-Fi"' (s')|^ in radiative return experiments (y^ = 
Mj.) for different angular cuts. Here the scenario is shown 
where the model error is estimated by a comparison of FSR 
from scalar particles with FSR from fermions of the same 
charge and mass [see (|3.33j)_] . The cut scenarios (a), (b) and 
(c) are the same as in Fig. m. 



where 



point,f/ N 



a^P^ cos ( s + imi 1 



4m^ 



- COS' 

3 







M 



(3.32) 

0^ being the minimal angle between the pion momenta 
and the beam axis allowed b y the given cuts. In analogy 
to the case without cuts in ( 3.2g ) we then may define a 
model error as 



JO),/ 



^"cutis') = 



(3.33) 



Fig. ^ shows the FSR contribution and Fig. the esti- 
mated model error for different kinematic cuts. It is inter- 
esting to note that the FSR contribution can be clearly 
reduced by the chosen kinematic cuts. This is especially 
obvious in the soft photon region. The suppression of the 
model error, being estimated by the scalar vs. fermionic 
scenario, by the considered cuts, however, is not obvious. 
As a matter of fact the cross sections drop out very quickly 
for low \/s' and vanish (up to higher order effects) below 
^/s' ~ 0.5 GeV (see Fig. 0). 

Finally we would like to comment on the estimation of 
FSR from a measurement of the pion forward-backward 
asymmetry |15 , defined as 

\ds') \ds') 



Afb(s') 



(da_ 

\ds' 



, cos e„<o 



' da 
<. ds' 



da 

^ \d7J 
.COS e„<o 

'int 



, cos 9^>0 



(da\ 
\ds') 



cos e„<o 



da 



cos 6I^<0 
'fin 



(3.34) 



Fig. 11. Pion pair invariant mass distribution in radiative 
return experiments (-^s = M<f) for the three discussed cut 



scenarios. 
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Fig. 12. Forward-backward asymmetry in per mill of nega- 
tively charged pions in radiative return experiments. is here 
the angle between the momentum of the incoming e~ and the 
outgoing n~ . 

Here 0^ is the angle between the momentum of the incom- 
ing and the outgoing 7r~. One could think that if Aps 
is small also the FSR contri bution must be suppressed, 
because, as can be seen from ( |3.34| ), the IFS contribution 
to the spectral function [da / ds')i-a± determines the size of 
ApB- The smallness of ApB: enhanced by kinematic cuts, 
could then be taken as a measure for the suppression of 
FSR in experiments. In Fig. ^v4i?B is shown for the three 
different cut scenarios. If we compare Fig. ^ with Figs. |^, 
^, H and |l^ we observe that we would get a different es- 
timation of model dependence from ApB than from our 
previous analysis of FSR. For example, for the case of no 
cuts we get an increasing model dependence from our in- 
vestigation of FSR if s' approaches threshold while Aps 
decreases for s' 4m?.. 
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Obviously, ApB is affected by FSR and hence can only 
be predicted by assuming a model like sQED. A compar- 
ison of such a prediction with experimental data is able 
to shed light on the validity of such a model. E.g., if we 
observe a good agreement of the measured ApB with the 
sQED prediction we could expect that likely also the FSR 
prediction by sQED could be a good approximation. 

Of course is not a direct measure of FSR and hence of 
the model dependence related to it. In fact the ratio of the 
FSR contribution and the IFS contribution is a strongly 
varying function of s' (see Fig. While in the region 
of the p resonance this ratio is relatively large it becomes 
small for low s'. This is true no matter whether cuts are 
applied or not. 



o 



V 



T3k3 



-(a): 


1 1 1 1 1 1 
no cuts 




: (b): 


7<'S0^S2O<', acse, / 




-(c): 




\(a) L 




^\ \ , 
\ ^, ^ 

(bV \ 

, , , , , w 


Ac) ~: 

■ 



0.4 0.6 0.8 1 

Vs' [GeV] 



Fig. 13. Ratio of the IFS contribution to final state contribu- 
tion for cos Ott < 0. Here 6tt is the angle between the momentum 
of the incoming e~ and the outgoing vr". 



error. The main reason is that without relying on ad hoc 
models like sQED we are not able to disentangle the to be 
extracted QCD quantity (pion form factor) from the real 
and virtual QED corrections. 

For radiative return measurements a direct extraction of 
a'^''^ (s) is not possible. In fact we can neither obtain a'^''^ (s) 
nor (t'*''(s) at 0{a) FSR precision without resorting to a 
model. Furthermore we have to deal with an 0(1) FSR 
background which is under control only if one of the fol- 
lowing criteria apply: i) it can be subtracted by using fac- 
torization in the soft photon region (s' ~ s), ii) it is sup- 
pressed by kinematic cuts (where it is possible) or iii) it 
is negligible in regions where a^^^{s') 3> o'^*'-'(s). We have 
shown that at <j> factories like DA<?NE we can control this 
FSR background only above Vs' ~ 500 MeV. 

We also had a look at the pion forward-backward asym- 
metry which is a model dependent quantity and thus is 
able to test model predictions against reality. Suppose 
that data would agree well with the sQED prediction this 
could be an indication that sQED also is able to describe 
FSR to some extent. However, it is not possible to esti- 
mate the FSR contribution in a straightforward way from 
a measurement oi ApB- The ratio of the FSR and the IFS 
correction is a strongly varying function of Vs'. Aps be- 
comes small for low •\/s' although the final state correction 
and hence the related model dependence is large. 

For radiative return measurements of \Fi^\s)\'^ at B- 
factories the 0(1) FSR background term is practically ab- 
sent. The model uncertainty due to FSR for the extraction 
of the undressed pion form factor in this case is determined 
by higher order effects. At BABAR the smallness of the 
observed cross section is compensated by a high luminos- 
ity. Initial state corrections here have to be known to a 
high precision. In particular for y/s' < 1 GeV the domi- 
nant pion pair production channel is e'^e~ 7r^7r~e"'"e~. 



4 Conclusions 



The importance of FSR for the extraction of the pion form 
factor from experimental data and for determining a^'^'^ 
has been discussed both for scan and radiative return ex- 
periments. 

We have shown that, by a photon-inclusive measurement 
and just subtracting ISR, a direct extraction of the FSR- 
inclusive cross section cr^'^^ in scan experiments is possible 
with excellent accuracy. In this case the model error is due 
to the breaking of ISR(8)FSR factorization by kinematic 
cuts and is of the order of a few per mill for the discussed 
C-symmetric angular cuts. 

On the other hand for exclusive measurements it is much 
more difficult to give a reliable estimation of the model 



In conclusion, to measure the FSR-inclusive pion form fac- 
tor \Fif\s)\'^ precisely and in a model independent man- 
ner at low energies, precise data from photon-inclusive 
scan experiments will be indispensable. We hope that such 
an experiment will be possible at DA<?NE at a later stage. 
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A Factorization of ISR and FSR 



Consider the process "e 



7 



X + photons" where 



X is an arbitrary non-photonic final state and the photons 
can either be emitted from the initial state (IS) or the final 
state (FS). We are interested in an expression for the in- 
clusive total cross section ct'"'^' = a{e^e~ photons) 
in terms of the FSR-inclusive cross section that is dressed 
by all real and virtual FS photonic corrections cr*^^^ and a 
universal IS radiator function pf^f corresponding to all IS 
real and virtual corrections. As will be shown in the fol- 
lowing such a factorization is in fact possible up to 0{a^) 
IFS real and virtual QED corrections. 

Let us consider first the process e+e^ — > 7* X -I- r 7, 
where r is a given number of real photons which can be 
emitted either from the IS or the FS. The amplitude A4^'^'> 
corresponding to this process can be written as the sum of 



all sub-amplitudes M 



{Vi,Vf) 



corresponding to real and 



Vi virtual photons attached to the IS e^e~ pair, rf real 
and Vf virtual photons attached to the final state X, and 
Vint additional virtual photons connecting the IS and the 
FS. For given r we have the condition Vi + rf = r. In 
the following we will neglect box- like diagrams, thus put 
Vint — 0. Hence we only keep the pure IS and FS vir- 
tual corrections (we will come back to the IFS contribu- 
tions later). The IR divergences of the virtual corrections 
are assumed to be regularized by a small photon mass. 
Without such an IR-regulator Ai'-^'' would not be defined. 
Obviously A^^*"^ by itself does not correspond to a phys- 
ical observable. However, at the end, after summation of 
all the contributions corresponding to all sub-amplitudes, 
we will obtain IR-finite, physical quantities and the IR 
regulator can be removed. The amplitude corresponding 
to the emission of r real photons can now be written as 



ri+r f—r 



00 



n+r f—r 



, (A.l) 



fcr'-*} are the IS real photon mo- 



where = {fcj'^ 

menta, {k'^^'^} = {k[^ ' . . .kVj'} the FS real photon mo- 
menta, {k^-^^} = {k''-^^ . . . fcix"*} the momenta of the non- 
photonic FS particles, q ~ pi + P2 the sum of the in- 
coming e"*"e~ momenta and finally qy ~ q — k^^^ = 
^ fc'"^^ -|- ^ fc*^-^-* the momentum of the virtua l ph oton 
■y*{qv) connecting the IS and the FS. In Eq. (A.l) we 



have written the sub-amplitudes M 



as contractions 



of rank-1 tensors A 



in 



containing the IS real and virtual 



corrections with the rank-1 tensors B, s containing the 
FS real and virtual corrections and the photon propagator 
related to "y*{qv)- This factorization of the amplitude is 



only possible because we have neglected the virtual IFS 
contributions. Defining 

00 

A(;'Hq,qv,{k^n) - E 4'^""'na,'?y,{fc«}) , 
B[;^)(qv,{fcW},{fc(^)}) = 



thus summing over all virtual IS and FS virtual correc- 
tions, we can write Ai^-^'^ simply as 



Y^B^^^^(qy,k^^\k^f^) 



(A.3) 



ri+rf—r 

Let us remember that Ai'-^^ contains all IS and FS vir- 
tual corrections to the given channel with r real (IS or 
FS) photons. Squaring the amplitude, averaging over the 
incoming e~^e~ spins and summing over the spins of the 
ISR photons S(i), the spins of the FSR photons S(j^) and 
the spins of the non-photonic particles S(x) yields 

'■^^''■'l^ i E EE E + IFS terms 

(A.4) 



= ^ |7W ('••■'•/) |2 + IFS terms 



ri,rf 

where "IFS terms" corresponds to real photon IFS contri- 
butions. Since we neglected already the virtual IFS con- 
tributions we have to do the same for the real IFS contri- 
butions. Otherwise we would have no cancellation of the 
IR divergences if later, after summing up all contributions 
r = ... 00, we want to remove the IR regulator. Neglect- 
ing the IFS contributions we can now express the squared 
amplitude |A4('''p as a sum of the incoherent contribu- 
tions |A^(''>'''/) p which can be written as a contraction of 
an IS and a FS tensor, respectively: 



\Mir.,rf)\2 ^ El:j\q,qv,{k^'^)F^:^^{qv,{k^''^}Ak^^'>}) , 

(A.5) 

with 

Ej[j\q,qv,{k^^^)= (A.6) 
i Al^'Hl, qv,{k^'^}) Ai^'^*{q, 9v, {fc«}) , 



(A.7) 



5:^B[:^,(gv^,{fcW},{fc(/)}) 
xB(-)(gv,{fc(^n,{fc'^^}) 
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Note that E^^[v\q, qv, {fc'*-*}) now only contains the IS real 
and virtual corrections while Fi^'^ . only contains the FS 
real and virtual corrections. Thus, already at this stage 
we obtain complete factorization of the IS and the FS 
corrections. Phase space integration over the r.i + rf + nx 
particle final state yields the total cross section related 
to the emission of r.i IS and ?"/ FS photons [(3e = (1 — 



^ri,r f 



is) 



1 



n rj nx 



n n n ^-^*p«a ^ dupsi/' dups. 



a=l b=l c=l 



{2.f s^'Kq-Y.^^^-Y.K'-Y.^i''') 



a=l 



b=l 



c=l 



with 



|A^('--'-/)|2 , 



(A., 



d^ki 



X) 



{2nf2E, 



(X)- 



(A.9) 



For the real photons the same IR regulator (photon mass) 
has to be used as for the virtual photons. Hence we are 
all the time deali ng with IR-regularized expressions. Al- 
though Eq. (A. 8) includes now virtual and real photon 
IS and FS corrections we cannot remove the IR regulator 
since we included only r real photons but virtual correc- 
tions to all orders, an^rjis) is obviously not a physical 
quantity by itself. Inserting the following identities (with 
sv = Qv), 



d'qv 5W(,^-^fc(/)-^fcW) 
b=l c=l 



and 

1 = y dsv S{sv - Qv), 
into Eq. (^J) yields 

<^T.,rf{s) = J dsv J n dLipsl^^ 

X 5(4)(5_g^-^fcW) 
a=l 

X 4;-)(g,gv,{fc^'^}).^C';)(9^) 

with the integrated FSR tensor 

„ »•/ nx 

^i^M^ = / n n d^'p'i^^ dupsi""^ 

■J U 1 1 



d^ 



2ql 



(A.ll) 



6=1 c=l 



' f nx 
l.(f) 



6=1 c=l 

F(';';)(g^,{fcW},{fc(/)}). (A.12) 



From Lorentz-covariance follows that the integrated FSR 
tensor can be written as a linear combination of the two 
linear independent tensors g^'^ and qyflv'- 

^A'^^fHsv)[g^''-'^^ (A.13) 



-tr 



Sv 



For the second equality in ( A.13 ) gauge invariance has 
been used, implying the Ward identity qv ^i^^^^-^iciv) = 0. 
Note that for the special case of no ISR photon emission 
[ri — 0) the related total cross section corresponding to 
the emission of r/ = r FSR photons can be written as 

-^.M-^^Ei^Kq)^^:,^{q), (A.14) 
with the lowest order IS tensor 

= [pl^.P2.+P2^.Plu - ^ g^") ■ (A. 15) 

Taking int o acco unt the Ward identity q^E^u) (q) — and 
using Eq. ( |A.13| ) we can write 



CTo,r/(s) = 



1 



2s/3e(s 
e'^ s + 2TOg 

2sPe{s) 



tr[4°)^(g)] ^tr [^(';^),(g) 



tr 



(A.16) 



For the general case including ISR photons the following 
Ward identity holds: 



<4;')(9,gy,{fc«}) = 



(A.17) 



Inserting the expression for T'^'^{qv) in ( A.13 ) into ( A.ll ) 
we can write the cross section corresponding to the emis- 
sion of ri ISR and rf FSR photons as 



(A.IO) crr.,r^(s) 



2sPe 



dsy — tr 
3 



n d^'P^a^ 



a=i ---y 

tr[M'--)''(g,gv,{fcW})] 



2q 



a=l 



(A.18) 



(note that tr ^(^j,yi^{qv) = 3A^^f\sv) is only a function 



of Sv)- Hence, using Eq. ( |A.16| ), we finally arrive at 



(^r.^rfis) = I dsv Cro,rf{sv) Pini\s,Sv) , (A. 19) 



with 



PkAis,sv) = - ^ 



1 



s + 2m^ 



a=l 

S^'Hq-qv-Yki:^) 

a=l 

t4El:'>{q,qv,{k^'^)] 



2g?. 



(A.20) 
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[Note that for the case of no ISR photons the above equa- 
tion directly giv es p^'f l (s, sy) = d{s—sv) which just projects 
out cro,r/ (s) in ( [A.19| )] . At this point we can ask the ques- 
tion what will be the expression for the inclusive total 
cross section (Tinci(s) with any number of real and virtual 
IS and FS photons. Neglecting the IFS contributions we 
can immediately write CTinci(s) as the incoherent sum of 
all ISR«)FSR contributions: 



Also the inclusive IS radiator we write as an expansion: 



oo oo oo 



Cilices) — ^ <^r,,rf{ 
Ti ,rf — 



(A.21) 

dsv a^^Hsv) sv) + 0(a2)iFs, 



with 



a^'^\sv) = (yG,ri{sv). fastis, sv) = Pini 

rj=0 ri=0 

(A.22) 



So in (A.21) we finally expressed the completely inclusive 
total cross section crinci(s) — crie^ e~ — > A" -I- photons) in a 
factorized form. Note that a'^'^^^sy) now contains the real 
and virtual FS corrections to all orders and p5"?'(s,sy) 
contains the real and virtual IS corrections to all orders. 
a'^^^i^s) and pjnf'(s, sy) are separately IR finite and the IR 
regulator can therefore be removed. This of course has to 
be the case since crinci(s) is [up to 0{p?) IFS effects] a 
physical observable. 

Using the above f ormul as it is now straightforward to ex- 
press iJinci in Eq. (A.21) as a perturbation series in a. For 
simplicity we will show explicitly the expansion only to 
0(a). For this we write the FSR-inclusive cross section as 
the expansion 



oo oo oo 



^^''Hs) = E E E ^(r,,.,..'^){s) , (A.23) 

r/=0 ''/=0 1)^=0 



with 



^ ^ '7 TLX 

'T(.„.„„^)(s) = tt^qt-: / n n <^^^vsS!^ dLipsi''') 

rf nx 

{2nr 5W(g - Y: ki^' - E k'c''') M;.,,.„„^),(A.24) 



and 



i>^i(.....;.) = -y(-^ + 2™DEE 



i3(,^,„^)^(g,{fcW},{fc(/)}) 



(A.25) 



rini 



(^'^w-EEE 

P{ri,Vi,v'-) {s,sv) (A.26) 

ri=0 Vi=0 v'.=0 



with 

P{ri,Vi,v'.)is, Sv) 



a=l 



A*;^„,)(q,gy,{fc''^}). (A.27) 
The inclusive cross section can then be written as 



CTinci(s)=y dsv<^ [ct(o,o,o)(sv) + o-(i^o,o)(sv) 

+ CT(o4,o)(sv) + o-{o,o,i){sv)] 

X [P{o,o,o)is,sv) + P(l,ofi)is,Sv) 

+P{o,ifi){s,sv) +P(o,o,i)(s,sy)]| +0(a^) 

= a(^)(s) [l + 6ini{s,A)] 

dsv Cr^^^(sy) Pini(s, Sy) + 0{a^), 

(A.28) 

with 

(J^'^Hsv) = o-(o,o,o)(sy) + cr(i_o,o)(sv) 

+(7(04,0) (sy) + CT(o,oa)(sy) + 0{a'^), 
Pini{s,Sv) = P(i,o,o)(s, Sy), 

Sini{s,A)^ / dsy [p(o,i,o)(s,sy) + P(o,o,i)('S,sy)] 



dsv P(i,o,o)(s,sy)- 



(A.29) 



Here A is the soft photon cut off, (7(o,o,o)(s) = "''■^Hs), ^^d 

P(0.0,0)(s, Sv) = S{s - Sy). 

The above derivation had many steps, however, it should 
be stressed that the most important one to get ISR(8)FSR 
factorization was th e use of Lorentz covariance and gauge 
invariance in (A. 13). If kinematical cuts are ap plied then 
ISR(8)FSR factorization breaks down because ( A. 13 ) will 
not be valid any more. Perturbatively this occurs already 
at 0{a), as we show in Appendix B. 



B Scalar and fermionic final state corrections 
to TT+TT production 



Taking final state corrections up to 0{a) into account the 
observed pion pair invariant mass distribution can be writ- 
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ten as the sum of an ISR and an FSR contribution: 

(B.l) 

For the analytic expression of the initial state radiator 
function pi„i(s, s'), including radiative corrections up to 
leading log 0{a^) and leading initial state e+e~ pair pro- 
duction contributions, we refe r to (17) in Integrating 
the spectral function in (B.l) over s' yields the observed 
total cross section (Tobs(s)- 

The 0(a) final state radiator functions, corresponding to 
hard photon radiation from scalar particles /5fin(s, s') and 
from fermionic particles Pf^^{s,s'), read {z = s'/s) 



Pfin(s, s') 



(B.2) 



X i?7r(s, s') [1 — Z 



<5(l-z) + 



(B.3) 



respectively, with the corresponding hard and soft photon 
functions 



dfi„(s,s)- ^ (1 Z)^3^^^ , 



i3^(s, s') 



a (1 — z) s P-k{s') 

( l + fi.{s') \ 

n s/3^(s) 



2/3. (s') 



■log 



I -Ms') 



(B.4) 



(B.5) 



(B.6) 



2 - 2 log 

2/3. 
log(/3.) 

2Li2 
2 o' 



log 



1 + 

1-/3. 



log 



log ( ^^f'' ] log 



2/3. 



1 + /3. 



2/3. 
2Li2 



1 + /3. 
2 

1-/3. 
2/3. 



1-/3. 
2/3. 



a 1 



7r2/3. 



1 + /3. 
1-/3. 



(B.7) 
(B.8) 



Neglecting soft photon exponentiation, we can write the 
observed total cross section as the sum of a soft photon 
contribution and a hard photon contribution: 



^o£(^)='^^"'(^) l + <5i„i(.,^) + <^(.,yl) 



(/)/ 



s~2-/sA 



4mJ 

a^°\s) 



ds' p,^,{s,s')a^'H-s') 



s-2\/syl 
4m?. 



(B.9) 



where 



Pfin(s, s') 

Sfin{s,A) = log 



'~Sg^is,s') + ^^^ 

~^"nAs,s') + ^j^ 



(B.IO) 
(B.ll) 

fin (S):(B.12) 



6Us,A) = log (^^j B.(.' = .) +5^+/(5).(B.13) 

Taking now only the leading order contribution to the 
total cross section and the final state corrections into ac- 
count leads to the 0{a) FSR~inclusive cross section 



= a(")(,s). 












ds' 










1 + 


IT 



(B.14) 



The analytic 0{a) expression for the integrated final state 
correction factors read 



vis) 



2Li2 - 



1+/3, 



-3 log (^) log (i±f-) - 21og(/3^) log (i 
-31og(^) -41og(/3.) 



i±A, 



for a scalar particle final state and 



(B.15) 



¥{s) = vis) + --k, 



s+2ml 2s/3. 



■^(5.-2™^^ 



-s/3.(s + 2m2) 



(.2 + 2™4)iog(i±|^) 
i_loe/'i±^^ 



1 



for a fermionic final state. 



(B.16) 



Finally some comments on applying kinematic cuts, lead- 
ing to a breaking of ISR(8)FSR factorization (see Appendix A). 
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Treating FSR by sQED we can write the spectral function 
with C-symmetric kinematic cuts as 



da 



(B.17) 



fda_ 



ini,cut 
point 



+ 0(^2), 



where {da / ds')^°™^^^. and {da / ds' are the correspond- 
ing IS and FS spectral function for point-like, scalar parti- 
cles. Note that within the considered sQED model we treat 
the non-perturbative QCD effects (pion form factor) in a 
factorized way which means that the FSR corrections are 
treated within pure sQED and therefore do not affect the 
pion form factor. By some straightforward manipulat ions 
we can write the total photon-inclusive cr oss se ction (2.3) 
for the given cuts in a similar form as in (A. 28). 



since in cr^^|(sy) sv corresponds to the invariant mass 
square of the pions including FSR, while in pl^l{s, sy) sy 
corresponds to the invariant mass square s' of the pions 
excluding FSR. 

The cut quantities are defined as follows: 

a(°i(s) = |^^i°)(s)|2a°r"'(s) (B.18) 

with 

a°ur'"(^) = - "'^%7^- (i - \ cos^ 

(B.I9) 

Here is the minimal angle between the pion momenta 
and the beam axis allowed by the given cuts. Similarly, we 
obtain 

a'^^^,{s)^\Fi^\s)\'alr-\s) (B.20) 



Let us note that in ( ) is a formal integration variable 



and 



f^{s,^) 



da\ 



point 



1 f da\ 



point 



0, point/ \ I Jt-' / 
^(cut) (■5) VCfS /fi„,(cut) 



(B.21) 
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